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ABSTRACT 

We present the results of an infrared spectroscopic survey of 23 late-M dwarfs with the NIRSPEC 
echelle spectrometer on the Keck II telescope. Using telluric lines for wavelength calibration, we are 
able to achieve measurement precisions of down to 45 m/s for our late-M dwarfs over a one to four 
year-long baseline. Our sample contains two stars with RV variations of >1000 m/s. While we require 
more measurements to determine whether these RV variations are due to unseen planetary or stellar 
companions or are the result of starspots known to plague the surface of M dwarfs, we can place upper 
limits of <40 Mjsini on the masses of any companions around those two M dwarfs with RV variations 
of <160 m/s at orbital periods of 10-100 days. We have also measured the rotational velocities for all 
the stars in our late-M dwarf sample and offer our multi-order, high-resolution spectra over 2.0 to 2.4 
lim to the atmospheric modeling community to better understand the atmospheres of late-M dwarfs. 
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1. INTRODUCTION 

With at least 70% of all stars within 10 pc identified as 
M dwarfs (Henry et al. 2006), these low-mass stars are 
becoming high priority targets for extrasolar planet sur- 
veys. With measurement precisions consistently at the 1 
m/s level, optical radial velocity (RV) surveys are capa- 
ble of detecting super-Earth mass planets around some 
of our nearest M dwarf neighbors. A classic example of 
the success of optical RV surveys around M dwarfs be- 
ing the GJ 581 system, which harbors up to six planets 
including a few with masses less than 10 Mjsini (Vogt 
et al. 2010; Mayor et al. 2008). The lightest known 
star with a planet is the M4V dwarf, GJ 1214, which 
has a 6.6 planet. This companion was discovered by 
the Mearth transit program (Nutzman & Charbonneau 
2008) and confirmed through optical RV measurements 
with a velocity amplitude of 12 m/s (Charbonneau et 
al. 2009). Despite M dwarfs comprising the majority of 
the super-Earth mass planets, they have been neglected 
in many of the original large scale RV programs because 
they are so optically faint. Out of the 185 M dwarfs 
within 10 pc, 35% of them are fainter than V=12 and 
16% of them are fainter than V=14 (T. Henry private 
comm.) A few surveys of early- to mid-M dwarfs with 
optical echelle spectrometers have achieved RV precisions 
of down to ~3 m/s and placed limits on the population 
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for periods of 10 to 100 days (Bonfils et al. 2012; John- 
son et al. 2010; Endl et al. 2006; Zechmeister et al. 
2006). These surveys, however, did not include any M 
dwarfs later than M6. One solution to the faintness is- 
sue inherent to optical observations of late-M dwarfs is 
to move the observations to the near-infrared part of the 
spectrum where M dwarfs are considerably brighter than 
solar-type stars. In addition, the contrast between the 
photosphere of the star and starspots is reduced at near- 
infrared wavelengths. These co-rotating features can be a 
prominent source of measurement noise when looking for 
the periodic variability produced by an orbiting planet. 
This is especially true for both M dwarfs and young stars 
which can have both flares and large starspots (Setiawan 
et al. 2008; Huelamo et al. 2008; Prato et al. 2008; 
Bailey et al. 2011). 

There are now a variety of infrared echelle spectrom- 
eters commissioned on 3-10 meter telescopes including 
NIRSPEC on Keck II (McLean et al. 1998), Phoenix on 
Gemini South (Hinkle et al. 2003), CRIRES on the VLT 
(Kaufl et al. 2004) and CSHELL on the IRTF (Greene et 
al. 1993). All of these spectrometers are capable of pro- 
ducing spectral resolutions of R>20000 within the 1.0 
to 2.5 micron wavelength region. In the near-infrared, 
the Earth's atmospheric spectrum is riddled with absorp- 
tion features from water, carbon dioxide and methane. 
This allows for a natural spectral reference frame from 
which to measure the change in the position of the stellar 
absorption lines (Griffin & Griffin 1973). Using telluric 
features, the CRIRES spectrometer achieved RV preci- 
sions of 5-20 m/s for bright M giant stars (Figueira et 
al. 2010; Seifahrt & Kaufl 2008). Other radial velocity 
precisions achieved with infrared echelle spectrometers 
using telluric calibration include 50-300 m/s for a sam- 
ple of nearby, slowly rotating early-M dwarfs (Bailey et 
al. 2012; Rodler et al. 2011; Crockett et. al. 2011) 
and 200 m/s for a sample of brown dwarfs (Blake et al. 
2010; Bailey et al. 2012). A precision of 50 m/s is suf- 
ficient to detect >2 Mj mass planets with periods of 10 
days around mid-M dwarfs. These infrared RV precisions 
can be compared to the ~ 1 m/s precision consistently 
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demonstrated around mid-M type stars at optical wave- 
lengths (Bonfils et al. 2011). 

Finding planets around M dwarfs is important because 
the two competing planet formation paradigms, core ac- 
cretion and gravitational instability, predict very differ- 
ent planet frequencies for the smallest stars: core accre- 
tion models predict relatively few planets (1%; Kennedy 
& Kenyon 2008) while gravitational instability models 
(e.g. Boss et al. 2006) predict frequencies similar to 
those of solar-type stars, perhaps 10% or more. Because 
M dwarfs are, by far, the most abundant spectral type in 
the stellar neighborhood, a thorough census of the plan- 
etary population around these stars is advantageous for 
potential follow-up studies. 

Here, we present the results of an infrared spectro- 
scopic survey of a sample of late-M dwarfs with the NIR- 
SPEC spectrometer on Keck II. In Section 2, we review 
the observing program used to collect the data as well 
as the data reduction and method used to measure the 
radial velocities of the spectra; in Section 3, we com- 
ment on the properties of the RV measurements for those 
stars with multiple measurements and stellar rotational 
velocities for the entire sample. In Section 4, we discuss 
the implications of these measurements and the future of 
planet searches with the NIRSPEC instrument. 

2. OBSERVATIONS AND ANALYSIS 

We have observed a sample of 23 late-M dwarfs over a 
spectral type range of M6.0 to M9.0 (see Table [1]) with 
the NIRSPEC spectrometer on the W. H. Keck II tele- 
scope (McLean et al. 1998). To avoid close binaries that 
would impede RV planet detections, these stars were cho- 
sen from high spatial resolution, adaptive optics surveys 
of nearby late-M dwarfs (Close et al. 2003; Seigler et 
al. 2003). Using the photometric distance relationship 
from Henry et al. (2006), we estimate the range in dis- 
tances for this sample of M dwarfs to be 4 to 50 parsecs 
based on their V-K colors and spectral types. Four of 
the M dwarfs have photometric distances placing them 
within 10 parsecs of the Sun. Seven of the M dwarfs 
were observed three or more times over a four year time 
span from 2006 June to 2010 July. GJ 752 A and B 
were added to our sample in mid-2009 due to the an- 
nouncement of an astrometric planet discovered around 
GJ 752 B (Pravdo & Shaklan 2009). NIRSPEC data 
for these two stars was collected by multiple observers 
during unrelated NIRSPEC runs with the same instru- 
mental setup. Additional K and early-M dwarfs (GJ 628, 
GJ 725 A, GJ 725 B and GJ 752 A, see Table p were 
observed as RV standards and to determine the infrared 
RV precision as a function of spectral type. Finally, a 
set of five A type stars were observed for determining 
the wavelengths solution of the instrument. 

We collected near-infrared spectra using the 
NIRSPEC-7 (2.22 ^im) passband in echelle mode 
with the 3 pixel slit (0.432") in combination with an 
echelle angle of 62.65 and grating angle of 35.50 degrees. 
This configuration results in a wavelength coverage of 
1.99 to 2.39 fim and a spectral resolution of R~24000 
(see Figure [j). Each epoch of observation included 
nodding the star along the slit to allow for subtraction 
of two images to remove background emission prior 
to spectral extraction. The integration times for the 
targets ranged from a few minutes per nod for the A 



stars and 30 minutes per nod for the late-M science 
targets. From the echelle data, which covers seven 
full spectral orders, we determine our radial velocities 
from the 33rd order on the focal plane covering the 
spectral range of 2.28-2.32 /xm. This wavelength region 
is ideal for RV measurements because it is permeated by 
telluric absorption features as well as the CO bandhead 
prominent in M dwarf spectra. 

2.1. Spectral Extraction and Modeling 

The data reduction and spectral analysis are described 
in detail in Bailey et al. (2012). Here we highlight the 
primary steps in the data reduction as well as any differ- 
ences between our analysis and that presented in Bailey 
et al. (2012). After flat fielding and nod-pair sky sub- 
traction, the 1-D spectra are extracted from the NIR- 
SPEC images using the 'optimal extraction' method first 
proposed by Home (1986). This method constructs spec- 
tra by summing weighted pixel values across the spectral 
profile, where the weights are based on noise statistics. 
Although the technique works well at removing deviant 
pixel values caused by such things as bad pixels and cos- 
mic ray events for high signal-to-noise spectra, occasional 
features sometimes remained in the lower signal-to-noise 
spectra of the faintest stars observed. These features 
appear as large positive or negative values of 1-2 pixels 
width and likely result from a poorly determined spec- 
tral profile at that location, perhaps caused by a larger 
than normal number deviant pixel values; this appears 
to be common for older generation infrared arrays. In 
these cases, these features were identified by eye in the 
spectra and the values were replaced by the median value 
of neighboring values. For each pair of observations, the 
spectra at each of the two nod positions are analyzed 
separately, providing two RV measurements per epoch. 

Each observed 1-D NIRSPEC spectrum is modeled as 
the composite of the intrinsic, rotationally broadened 
spectrum of the star and the telluric spectrum convolved 
with the instrumental profile. The telluric spectrum is 
extrapolated from spectral observations of the Sun (Liv- 
ingston & Wallace 1991) while the intrinsic spectrum 
of each of the late-M dwarfs is extrapolated from the 
NextGen models, which account for the star's effective 
temperature and surface gravity (Hauschildt et al. 1999). 
An initial estimate of the wavelength calibration and 
spectral response function are determined from a set of 
A star calibrators observed during the same night as the 
science targets. A stars are ideal spectral calibrators be- 
cause the spectrum of the star is featureless in the wave- 
length region used to estimate the radial velocities. To 
determine the initial values of the wavelength solution 
and instrumental spectral profile, the model telluric spec- 
trum is fitted to the observed spectrum of the A star by 
minimizing the \ 2 value with the amoeba downhill sim- 
plex algorithm (Press et al. 1992). The fitting process 
involves four free parameters: three for the quadratic 
wavelength solution and one for the gaussian FWHM of 
the instrumental spectral profile. 

With initial guesses taken from the known properties 
of the spectrograph determined by previous NIRSPEC 
spectrograph programs, all four free parameters are al- 
lowed to vary until the best fit is achieved. The four pa- 
rameters characterizing the instrument are then used as 
initial guesses for finding the best fitting model spectrum 
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of the science targets. There are ten free parameters 
involved when modeling the observed M dwarf spectra: 
three for the quadratic wavelength solution, and one each 
for the instrumental spectral profile, limb darkening, RV, 
vsini, and airmass. There are also two scaling factors for 
the relative depth of the telluric and synthetic M dwarf 
spectrum. When available we use published values of the 
vsini for the initial guess during the fitting process (see 
Table IU. 

To obtain the optimal fit to the observed M dwarf 
spectra, we vary the wavelength solution parameters, the 
FWHM of the instrumental spectral profile, and the ra- 
dial velocities and normalization parameters. The fitting 
is stopped when the x 2 value ceases to change between 
variable perturbations by more than 10~ 3 %. Figure [2] 
shows the best fit to one epoch of the spectra for the 
slowly rotating M3.5 and rapidly rotating M7.5 dwarfs, 
GJ 752 B and 2M1757+70. The telluric spectrum and 
the theoretical stellar spectrum rotationally broadened 
and matched to the published stellar spectral type are 
also plotted. 

In some cases, a disagreement in the estimated radial 
velocity determined for the A and B nod pairs (> 100 
m/s) during a single epoch requires the manual pertur- 
bation of the initial guesses to the fitting parameters in 
order to bring the two nod RV values into agreement. 
This additional step was necessary because the amoeba 
algorithm is sensitive to the initial guesses. We adopted 
a methodology that automatically finds the appropriate 
initial guesses but still requires a small amount of cus- 
tomized adjustment by hand to find the best fit to the 
data. For each epoch of observation we calculate the 
average of the radial velocities determined from all spec- 
tral nod pairs. As described in Bailey et al. (2012), the 
observed RV precision over multiple epochs is the combi- 
nation of the theoretical photo noise error, NIRSPEC's 
instrumental error and the error caused by intrinsic stel- 
lar activity. The theoretical photon noise error derived 
assuming a perfect instrument, no intrinsic stellar jitter, 
the signal-to-noise of the observed spectrum and vsini of 
the star (Butler et al. 1996). Our NIRSPEC observa- 
tions of the RV standards (GJ 628, GJ 725 A and GJ 
725 B) result in RV precision of 75-143 m/s for nearby, 
bright, slowly rotating mid-M stars (Bailey et al. 2012). 
Observations of these same RV standards with greater 
instrumental precision suggest that they have an intrin- 
sic stellar noise error of < 10 m/s. The uncertainties 
quoted for a single RV measurement (see Table [3]) arc 
then the quadrature sum of the photon noise error (varies 
by star), instrumental error (103 m/s, the average of the 
RMSs from the RV standards) and stellar jitter error (10 
m/s). In Table H] we also include the error associated 
with each data point (a x 2 =1 ) that would produce a x 2 
of one if there was no detectible planet around the star 
resulting in a flat RV curve. 

3. RESULTS 

The observed standard deviations of the RV measure- 
ments for those stars with at least three epochs of ob- 
servations are listed in Table |4] and range from 45 to 
2015 m/s for those seven stars in our primary sample 
of 23 late-type M dwarfs with three or more observa- 
tions. There are two stars (2M2313+21 and 2M0027+22) 
that have RV dispersions greater than 1000 m/s and two 



stars (2M0320+18 and 2M1757+70) with RV dispersions 
greater than 500 m/s. While the fact that one star has 
an RV dispersion of 45 m/s (2M0253+16) could repre- 
sent the optimal precision available with this observing 
method, this comes from only three individual measure- 
ments and should be considered tentative. 

The M6 star, 2M2313+21, has the largest scatter in its 
RV measurements (2015 m/s) and an intermediate rota- 
tional velocity (16 km/s). However, with only three RV 
observations, its difficult to conclude whether the large 
change in the velocity of the star from 2006 to 2010 is 
due to an unseen companion. The same issue is also true 
for 2M0320+18 which has an intermediate observed ro- 
tational velocity of 15 m/s and an RV dispersion of 546 
m/s and 2M0027+22 which has the highest observed ro- 
tational velocity (60 km/s) of all the stars in the sample 
and an RV dispersion of 1050 m/s. The combination of 
too few RV measurements, late-type spectral types and 
large (> 1 km/s) rotational velocities makes the clas- 
sification of any of these stars as companion candidate 
premature in the absence of additional RV measurements 
and a larger sample size. Additional RV measurements 
will allow us to discern between RV variability due to 
companions, starspots or the measurement floor as a 
function of stellar spectral type from the NIRSPEC tel- 
luric calibrations. 

The M7.5 dwarf, 2M1757+20, has six RV epochs mea- 
sured over the span of four years (see Table |4j . From 
these observations we see a standard deviation in the 
RV measurements of 573 m/s (see Figure |3|). Using 
the theoretical RV precision, number and cadence of the 
NIRSPEC observations, we utilize Monte Carlo simula- 
tions to place limits on the presence of planets around 
2M1757+20. Because this survey is most sensitive to 
short period companions, we set limits at representa- 
tive orbital periods of 10 and 100 days. A set of 10,000 
circular orbits are generated with random inclinations, 
phases, and companion masses ranging from ~1 to 80 
Jupiter masses. Circular orbits are a valid assumption 
for the 10 day orbits if we are assuming that the planets 
were subject to circularization via disk interactions (Lin 
et al. 1996). For the 100 day orbits, the addition of or- 
bits over multiple eccentricities effect our mass limits in 
a way that is smaller that our current measurement un- 
certainties. In combination with the stellar masses listed 
in Table [IJ these orbits are used to produce a set of RV 
measurements with the same temporal sampling as our 
NIRSPEC observations. The values of the model RV 
measurements are perturbed by an error value randomly 
chosen from a Gaussian distribution with a FWHM of 53 
m/s. This observational error corresponds to the median 
error due to the photon noise and instrument error. Con- 
servative detection limits are then set by determining the 
companion mass that induces an RV dispersion that is 2a 
greater than the photon error of the star depending on its 
rotational velocity. Using a maximum RV amplitude of 
573 m/s for 2M1757+20 and assuming a circular orbit, 
the companion responsible for this perturbation would 
have a mass of <16 Mjsini in a 10 day orbit and <40 
Mjsini in a 100 day orbit. These limits do not rule out 
that the companion could be a brown dwarf or very-low 
mass star. 
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3.1. GJ 752 B 

Astrometric measurements of the late-M star GJ 752B 
(a.k.a. VB 10) by Pravdo & Shaklan (2009) caused great 
excitement as they showed perturbations that were in- 
terpreted as an unseen planet with a mass of 6.4 Mj 
and period of 0.74 years. This result was touted as the 
first exoplanet found solely through astrometric measure- 
ments. However, the discovery was short-lived as sub- 
sequent radial velocity measurements quickly proved the 
planet to be non-existent (Bean ct al. 2009) and more re- 
cent astrometric observations provided further evidence 
that the original astrometric detection suffered from sys- 
tematic errors (Lazorenko et al. 2010). Our observations 
continue the tradition as we observed GJ 752 A and B 
eight times from June 4 2009 to October 31 2009 (see 
Table [3]). The dispersions of our nightly RV measure- 
ments are 143 m/s and 351 m/s for GJ 752 A and GJ 
752 B, respectively (see Figure [3]). The scatter in our 
NIRSPEC RV measurements are smaller than the mea- 
surements of Zapatero Osorio et al. (2009) which reached 
an amplitude of 1 km/s. Our RV dispersion for GJ 752 B 
is larger than the 10 m/s reached by Bean et al. (2009) 
and could be due the smaller spectral resolution of NIR- 
SPEC compared to CRIRES (R=24000 for NIRSPEC 
versus R=50000 for CRIRES). Monte Carlo simulations 
of our RV measurements place limits on the mass of a 
planet around GJ 752 B of 7 and 14 Mj for periods of 
10 and 100 days, respectively. 

3.2. Rotational Velocities of Late-M dwarfs 

Because the stellar rotational velocity is one of the free 
parameters used when modeling the observed spectrum, 
we are able to determine this physical parameter for all of 
the stars in our late M dwarf sample, which range from 6 
to 60 km/s (see Table 0J. Figure 2] plots the values of the 
rotational velocities estimated from our infrared spectra 
as a function of spectral type along with those values for 
a sample of early-L dwarfs also observed with NIRSPEC 
(Blake et al. 2010) as well as late-M dwarfs observed 
at optical wavelengths (Mohanty & Basri 2003; Reiners 
& Basri 2010). The uncertainty in the rotation velocity 
is determined by running the x 2 minimization program 
while only allowing the vsini to vary. Our stated uncer- 
tainties are then the standard deviation of all the result- 
ing rotational velocities that provide the best fit to the 
observed spectra. For those stars with only one epoch 
of observations (two spectra), the uncertainty stated is 
the half range of the vsini values. Our infrared vsini val- 
ues agree with those estimated from optical spectra and 
recently published J (1.25 /xm) NIRSPEC spectra (Rein- 
ers & Basri 2010; Deshpande et al. 2012, see Table g]). 
When we combine our vsini values with those from the 
literature (see Figure 2J , our values lie squarely within 
the spread of values from previous studies (Blake et al. 
2010; Mohanty & Basri 2003; Reiners & Basri 2010), 
which show a trend of increasing minimum rotational 
velocity with later spectral type. 

4. DISCUSSION AND CONCLUSIONS 

While we are not able to confirm the detection of 
any planetary or stellar mass companions with this pilot 
study of 23 late-M dwarfs with the NIRSPEC instrument 
on Keck, we have demonstrated that we are capable of 



producing RV dispersions of down to 50 m/s for late-M 
dwarfs with intermediate rotational velocities and telluric 
wavelength calibrations. Additional stars and RV mea- 
surements are needed to definitively conclude whether 
those stars with large RV dispersions have unseen com- 
panions or whether we are encountering star spots or 
the noise floor for these spectral types and instrument 
configuration. With radial velocity precisions down to 
~50-150 m/s for stars of masses of 0.08-0.15 M Q , these 
observations would be sensitive to planets with masses of 
<1 Mjsini at periods of 10-100 days. Considering there 
are a few thousand known M dwarfs within 25 parsecs of 
which a few hundred have V-Ks colors suggesting they 
are M5 spectral type and later (T. Henry private commu- 
nication), and that these stars represent a population of 
stars largely unexplored by most planet detection tech- 
niques, these observations open the door to potentially 
fruitful new discoveries. While the percentage of early- 
M dwarfs with gas giant planets is currently estimated 
at 3% (Johnson et al. 2010), we have tentative informa- 
tion that this percentage could be < 1% for late- type M 
dwarfs at the smallest separations of a<0.05 AU (Blake 
et al. 2010) and little information on these planetary 
populations at larger separations. With the masses of 
these stars approaching that of the most massive brown 
dwarfs, the planet population around them could have 
significant implications on the method of planet forma- 
tion that is dominant in this stellar mass regime. 

Saar et al. (1998) and others have used starspot mod- 
els and observations to show that the radial velocity jitter 
of a star is directly related its rotational velocity. There- 
fore, determining the rotational velocities for a large por- 
tion of all the of the nearby M dwarfs and determining 
whether there is any observed correlation between stellar 
rotational velocity and optical or infrared radial velocity 
precision would greatly benefit future planet search ef- 
forts. 

These high-resolution NIRSPEC spectra can be used 
for more than just radial velocity and vsini measure- 
ments. Our data represent a additional component of a 
small sample of K-band, echelle spectra of late-M dwarfs 
published to date (i.e. Blake et al. 2010). While our 
radial velocity program uses just the 2.28-2.31 /iin order, 
the NIRSPEC data cover a wavelength range of 1.98-2.39 
zxm (see Figure [T]). With the discovery of brown dwarfs 
a few decades ago, modeling the spectral features and, 
therefore, the complex atmospheres of late-M dwarfs has 
not been as thorough over all near-infrared bands when 
compared to what has been accomplished for solar-type 
stars and brown dwarfs. None of our spectral show signs 
of Br 7 emission (2.166 /zm) which could be used to aid 
in age determinations, however, fully reduced versions 
of our A star and M dwarf spectra are available to any 
who wish to utilize them for spectral synthesis and RV 
templates. 

With so few targets in this pilot study, it is difficult 
to determine whether there is any observational corre- 
lation between radial velocity precision and the stellar 
rotational velocity or spectral type of the star. How- 
ever, we have demonstrated that, with additional obser- 
vational epochs and a larger sample size, we could de- 
tect hot-Jupiter planets around late-M dwarf stars, which 
are often too faint to be observed efficiently with optical 
echelle instruments. Planned upgrades to the NIRSPEC 
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instrument , which include adding a methane gas absorp- 
tion cell to improve wavelength calibration, could lead 
to radial velocity precisions of ~20-30 m/s (Bean et al. 
2009). Such observations would be sensitive to Jupiter- 
mass planets in orbits of less than a year around late-M 
dwarfs. With an emphasis in detecting Earth- like plan- 
ets set forth by the 2010 Decadal Survey and in searching 
for planets around M dwarfs set forth by the 2008 Exo- 
planct Task Force, radial velocity surveys with infrared 
echelle spectrometers will have a significant impact on 
future exoplanet studies. 
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TABLE 1 
Sample of Late M dwarfs 



Target 


RA 


Dec 


SpTy 


K s 


2MUU27+22 


nn ot c tr n 

00 27 55.9 


i oo i n oo o 

+22 19 o2.8 


M8.0 


9.57 


o"\/Tm /in i nv 

2MU14U+27 


m ac\ no 

01 40 02.6 


i ot n 1 tr n c 

+27 01 50.6 


M8.5 


11.43 


2MU227-lb 


no OT -\ C\ A 

02 27 1U.4 


1 C O A A o n 

-lb 24 4o.U 


M9.0 


12.14 


oATnntro i 1 r* 

2MU25o+lb 


no rr o nn n 

02 53 00.9 


i 1 c tr o tro o 

+ 16 52 53.3 


M6.5 


7.59 


oA/rnntro i 07 

2MU25o+27 


no c o on o 

U2 5o 2U.o 


1 07 1 O O O O 

+27 lo 66.2 


M8.0 


11.48 


OA/Tnon^; oa 
2MUoUb-ob 


no r\c 1 i c 

Uo Ub 11. b 


na AT CO n 

-6b 4 ( 52.9 


M8.5 


10.63 


oa /f noon i i q 
2MUo2U+lo 


no on en 7 

Uo 2U 59. 7 


1 1 Q rr, A OO Q 

+ 18 54 2o.o 


M9.0 


10.64 


2MU /52+lb 


(17 to OO n 

()( o2 zo.y 


1 1 C 1 O 1 C T 

+ lb 12 15.7 


M / .0 


9.85 


oatooio 1 no 

2MUolo+2o 


no "i o co 1 

08 18 58.1 


i oo o o tr o o 

+2o oo 52.2 


M7.0 


11.15 


2M0853-03 


no cro or* o 

08 53 3d. 2 


no on oo 1 

-Uo 29 o2.1 


M9.0 


9.94 


ZM1U49+Z0 


in /in a i c: 
1U 49 41.0 


i O c QO co c 

+25 oo 5o.b 


M6.0 


11 A A 

11.44 


2M1124+13 


11 24 53.3 


+13 22 53.4 


M6.5 


10.07 


2M1546+37 


15 46 05.4 


+37 49 45.8 


M7.0 


11.41 


2M 1707+64 


17 07 18.3 


+64 39 33.1 


M9.0 


11.38 


2M1757+70 


17 57 15.4 


+70 42 01.2 


M7.5 


10.40 


2M1835+32 


18 35 37.9 


+32 59 54.6 


M8.5 


9.17 


GJ752 B a 


19 16 57.6 


+05 09 02.0 


M8.0 


8.77 


2M2052-23 


20 52 08.6 


-23 IS 09 6 


M6.5 


11.29 


2M2221+27 


22 21 54.4 


+27 29 07.0 


M6.0 


11.52 


2M2235+18 


22 35 49.1 


+18 40 29.9 


M7.0 


11.37 


2M2306-05 


23 06 29.3 


-05 02 28.6 


M7.5 


10.30 


2M2313+21 


O O "1 O A *7 A 

23 13 47.4 


+21 17 29.4 


M6.0 


10.44 


2M2349+12 


23 49 49.0 


+12 24 39.0 


M8.0 


11.56 


a aka VB 10 














TABLE 2 










Calibration Stars 






Target 


RA (2000) 


Dec (2000) 


SpTy 


K s 


A stars 


HR 104 


00 28 13.7 


+44 23 40.0 


A2 


9.99 


HR 8518 


22 21 39.4 


-01 23 14.4 


AO 


4.02 


HR 5511 


14 46 14.9 


+01 53 34.4 


AO 


3.65 


HR 7724 


20 14 16.6 


+15 11 51.4 


AO 


4.77 


HR 7336 


19 20 35.7 


-00 53 31.8 


B9 IV/V 


5.46 


RV Standards 








GJ 628 


16 30 18.1 


-12 39 45.3 


M3.5 


5.08 


GJ 725 A 


18 42 46.7 


+59 37 49.4 


M3.0 


4.43 


GJ 725 B 


18 42 46.9 


+59 37 36.7 


M3.5 


5.00 


GJ 752 A 


19 16 55.3 


+05 10 08.1 


M3.5 


4.67 
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TABLE 3 
Radial Velocity Measurments 



Star 


HJD a - 2400000 


Date [Y.M.D] 


RV [m/s] 


2M0027+22 


53928.576 


2006.07.12 


-9633±181 




55401.559 


2010.07.24 


-11428±226 




55402.555 


2010.07.25 


-11475±202 


2M0253+16 


53928.620 


2006.07.12 


68330±107 




55401.630 


2010.07.24 


68253±106 




004U2.DO 7 


om n o7 ok 
2U1U.U/.ZO 


CQOKnxi no 


2M0320+18 


53931.633 


2006.07.15 


45953±118 




55401.610 


2010.07.24 


45638±111 




55402. ozl 


o r\ i r\ rev o tr 

2010.07.25 


a a onn I i n a 

44o9U±104 


2M1707+57 


53929.415 


2006.07.13 


-10001±151 




54308.372 


2007.07.27 


-9990±185 




tr /i onn inn 

54309. 390 


onnT rev oo 

2007.07.28 


-11184±149 




54311.367 


2007.07.30 


-11042±150 




54312.428 


2007.07.31 


-11153±147 




55401.452 


2010.07.24 


-11010±181 


2M2306-05 


53929.583 


2006.07.13 


-52919±108 




55401.536 


2010.07.24 


-52600±106 




55402.581 


2010.07.25 


-52746±106 


2M2313+21 


53929.562 


2006.07.13 


-7391±116 




55401.582 


2010.07.24 


-10873±117 




55402.602 


2010.07.25 


-10889±118 


GJ 752b 


55016.448 


2009.07.04 


35614±104 




55024.496 


2009.07.12 


35644±103 




55041.259 


2009.07.29 


35738±102 




55110.281 


2009.10.06 


36483±102 




55115.297 


2009.10.11 


36161±102 




55116.315 


2009.10.12 


36169±104 




55135.202 


2009.10.31 


36491±115 




55401.510 


2010.07.24 


36081±105 



HJD - heliocentric Julian date. All measurements heliocentric radial velocities. 
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TABLE 4 

RV RMS and vsini Values for RV Standards and Science Targets 



Target 


N 


RV 


RV RMS 


<V i 


vsini - this work 


vsini - published" 






km/s 


m/s 


m/s 


km/s 


km/s 


f~< T £JOQ 


1U 


-Zl.l 


on 
oU 






i 
i 


(~* T 79^ A 





-u.o 


CO 

oy 






i 
i 


/~< t 7oc; T3 

t^rj (ZO D 


1 Q 


1 /i 
1.4 


1 






i 


CjrJ (52 A 


5 


36.2 


1 /I Q 

143 




... 


2 




8 


36.0 


351 


350 


8±7 


8 


2MUUz7+zz 


3 


-10.8 


1050 


820 


60±12 


56 


2M0140+27 


1 


10.4 






7±4 


<12 


zlYlUzz7-lo 


2 


49.0 






17±11 




2M0253+16 


3 


68.3 


45 


37 


10±4 




2M0253+27 


1 


47.8 






18±5 


17.5 


oi\,Tr\or\(^ or 

2MUol)o-oo 


2 


12.2 




.*** 


10±7 


■■■ 


zMUozU+lo 


3 


45.5 


04b 


430 


15±6 


15 




1 


-15.5 






6±4 




2M0818+2 3 


1 


33.5 






8±3 




oi\/Tr\ocro r\o 

2MUo5o-Uo 


1 


10.3 






20±7 






i 
i 


/ .O 






s-t- 
ozto 




2M1124+13 


1 


8.4 






7±4 




2M1546+37 


1 


-22.5 






10±7 


10 


2M1707+64 


2 


-10.1 






33±12 




2M1757+70 


6 


-12.3 


573 


500 


34±11 


33 


2M1835+32 


2 


4.8 






34±17 


44 


2M2052-23 


1 


35.8 






18±10 




2M2221+27 


1 


-14.8 






6±4 




2M2235+18 


1 


-14.9 






16±8 




2M2306-05 


3 


-52.8 


160 


130 


8±5 


6 


2M2313+21 


3 


-9.3 


2015 


1600 


16±4 




2M2349+12 


1 


-2.0 






12±4 


<12 



a Rcincrs & Basri 2009 (optical) and Dcshpandc ct al. 2012 (infrared). 
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2.04-2.07 
10 1- 2.06-2.13 
2.15-2.17 



2.22-2.25 
2.28-2.32 
2.36-2.39 



Fig. 1. — Plot of all orders of GJ 752 B observed (with wavelengths labeled in microns) with the NIRSPEC-7 filter. In addition to using 
the 2.28-2.32 (im region for our RV determinations we can also use this data to search for Br 7 emission and conduct spectral modeling to 
determine the properties of the atmosphere of these stars. These spectra still contain features due to telluric absorption. 
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atmospheric transmission spectrum 
theoretical spectrum 



observed spectrum and best fit model (SNR: 243) 



residuals (mean: 3.25%) 



f.,r-~ 
2.290 



2.295 2.300 2.305 

Wavelength (microns) 



atmospheric transmission spectrum 
theoretical spectrum 

observed spectrum and best fit model (SNR: 102) 
residuals (mean: 3.15%) 



2.295 2.300 2.305 

Wavelength (microns) 



Fig. 2. — Plot of the different components involved in modeling the spectrum of the slowly rotating (vsini = 8 km/s) M3.5V dwarf GJ 752 
B (left) and and the rapidly rotating (vsini = 61 km/s) M7.5V dwarf 2M1757+70 (right) for the purpose of determining the radial velocity 
of the star. They include the telluric spectrum (top) extracted from a high resolution solar spectrum, the spectrum of the target star as 
determined from NextGen models (middle), and the combination of these as well as an assumed vsini and instrument profile (bottom). 
The bottom dotted line is the difference between the best fitting model and the data. 
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500 1000 
JD - 2453522 [days] 



1000 : 
500 : 



2M 1757+70 



-500 \ 
-1000 | 



500 1000 
JD - 2453929 [days] 



1000 ; 
500 : 



GJ 752b 



-500 ; 
-1000 : 



500 1000 
JD- 2455016 [days] 



Fig. 3. — Top: Plot of the RV measurements for the standard star GJ 628. Middle and Bottom: Plots of the RV measurements for GJ 752 
B and 2M1757+70. The error bar on the top, left corner of each plot is scaled to the average of the errors of each of the RV measurements 
which are derived from the photon errors (50 - 100 m/s). 
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Fig. 4. — Plot of our estimated rotational velocities as a function of stellar spectral type (large diamonds) along with those determined 
for a sample of early L dwarfs also with NIRSPEC (Blake et al. 2010, plusscs) and additional late-M dwarfs with vsinis determined from 
optical echelle spectra (Mohanty & Basri 2003, Reiners & Basri 2010, triangles). Our data agrees with previous studies in that there 
appears to be an increase in the rotational velocity with later spectral type. 



